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Abstract 

We consider the problem of extracting a maximum-size reflected network in a 
linear program. This problem has been studied before and a state-of-the-art SGA 
heuristic with two variations have been proposed. 

In this paper we apply a new approach to evaluate the quality of SGA. In 
particular, we solve majority of the instances in the testbed to optimality using a 
new fixed-parameter algorithm, i.e., an algorithm whose runtime is polynomial in 
the input size but exponential in terms of an additional parameter associated with 
the given problem. 

This analysis allows us to conclude that the the existing SGA heuristic, in fact, 
produces solutions of a very high quality and often reaches the optimal objective 
values. However, SGA contain two components which leave some space for im- 
provement: building of a spanning tree and searching for an independent set in a 
graph. In the hope of obtaining even better heuristic, we tried to replace both of 
these components with some equivalent algorithms. 

We tried to use a fixed-parameter algorithm instead of a greedy one for search- 
ing of an independent set. But even the exact solution of this subproblem improved 
the whole heuristic insignificantly. Hence, the crucial part of SGA is building of a 
spanning tree. We tried three different algorithms, and it appears that the Depth- 
First search is clearly superior to the other ones in building of the spanning tree for 
SGA. 

Thereby, by application of fixed-parameter algorithms, we managed to check 
that the existing SGA heuristic is of a high quality and selected the component 
which required an improvement. This allowed us to intensify the research in a 
proper direction which yielded a superior variation of SGA. This variation signif- 
icantly improves the results of the basic SGA solving most of the instances in our 
experiments to optimality in a short time. 
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1 Introduction, terminology and notation 



Large-scale LP models which arise in applications usually have sparse coefficient ma- 
trices with special structure. If a special structure can be recognized, it can often be 
used to considerably speed up the process of solving the LP problem and/or to help in 
understanding the nature of the LP model. A well-known family of such special struc- 
tures is networks; a number of heuristics to extract (reflected) networks in LP problems 
have been developed and analyzed, see, e.g., [3 6,|5]|7][T2][T3][T8) (a formal definition 
of a reflected network is given below). From the computational point of view, it is 
worthwhile extracting a reflected network only if the LP problem under consideration 
contains a relatively large reflected network. 

We consider an LP problem in the standard form stated as 

Minimize {p T x; subject to Ax — b, x > 0}. 

LP problems have a number of equivalent, in a sense, forms that can be obtained from 
each other by various operations. Often scaling operations, that is multiplications of 
rows and columns of the matrix A of constraints by non-zero constants, are applied, 
see, e.g., J3][6][7][T2). In the sequel unless stated otherwise, we assume that certain 
scaling operations on A have been carried out and will not be applied again apart from 
row reflections defined below. A matrix B is a network (matrix) if B is a (0, ±\)-matrix 
(that is, entries of B belong to the set {1,0,-1}) and every column of B has at most 
one entry equal to 1 and at most one entry equal to —1. The operation of reflection of 
a row of a matrix B changes the signs of all non-zero entries of this row. A matrix B 
is a reflected network (matrix) if there is a sequence of row reflections that transforms 
B into a network matrix. The problem of detecting a maximum embedded reflected 
network (DMERN) is to find the maximum number of rows that form a submatrix B of 
A such that B is a reflected network. This number is denoted by v(A). The DMERN 
problem is known to be NP-hard [4). 

Giilpinar et al. (T3 1 showed that the maximum size of an embedded reflected net- 
work equals the maximum order of a balanced induced subgraph of a special signed 
graph associated with matrix A (for details, see Section|2]i. This result led Giilpinar et 
al. Ifl3l to a heuristic named SGA for detection of reflected networks. Computational 
experiments in lfT~3l with SGA and three other heuristics demonstrated that SGA and 
another heuristic, RSD, were of very similar quality and clearly outperformed the two 
other heuristics in this respect. However, SGA was about 20 times faster, on average, 
than RSD. Moreover, SGA has an important theoretical property that RSD does not 
have: SGA always solves the DMERN problem to optimality when the whole matrix A 
is a reflected network [ 13"]. Since SGA appeared to be the best choice for a heuristic for 
detection of reflected networks, Gutin and Zverovitch lTT4l investigated 'repetition' ver- 
sions of SGA and found out that three times repetition of SGA (SGA3) gives about 1 % 
improvement, while 80 times repetition of SGA (SGA80) leads to 2% improvement. 

In this paper we propose a more refined analysis of the SGA heuristic. By using 
a fixed-parameter algorithm, we managed to find the optimal solutions for majority of 
the instances in the testbed. This helped us to understand that SGA, in fact, obtains 
solutions of very high quality and sometimes even solves the instances to optimality. 
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This means that even a small improvement of SGA quality is a significant achievement. 
Having this result, we tried to improve SGA. 

SGA contains two components which leave some space for improvement. One is 
an independent set searching algorithm. In the original version, a greedy algorithm 
was used for this purpose. We replaced it with a fixed-parameter algorithm. Here we 
used the well-known fact that the complement of an independent set in a graph is a 
vertex cover. However, our experiments have shown that even the greedy algorithm 
usually finds the optimal or almost optimal solutions and, thus, this modification of 
SGA appears to be of little practical interest. This demonstrated that the independent 
set extracting heuristic need not be replaced by a more powerful heuristic or exact 
algorithm. 

Hence, the crucial component in SGA is the second component, i.e., building of a 
spanning tree. In the original version of SGA ( lfT3l[T4l ) a random procedure was used 
for this purpose. We tried to replace it with Breadth-First Search (BFS) and Depth- 
First search (DFS) algorithms. The experiments show that the final solution quality 
significantly depends on spanning tree and that using DFS is clearly superior to both 
BFS and Random Search (RS). 

Moreover, we observe that the choice of the algorithm for spanning tree compu- 
tation does not essentially influence the runtime of the heuristic under consideration. 
Thus we propose three new heuristics SGA(DFS), SGA3(DFS), and SGA80(DFS) that 
take about the same time as their respective counterparts SGA, SGA3, and SGA80 but 
are more precise. 

In order to evaluate the quality of the considered heuristics we compare their out- 
puts with optimal solutions of the considered instances. To solve the instances to op- 
timality, we design a fixed-parameter algorithm for the DMERN problem. Such algo- 
rithm is polynomial in the input size but exponential in terms of an additional parameter 
associated with the given problem. Problems that can be solved by fixed-parameter al- 
gorithms are called fixed-parameter tractable (FPT). A fixed-parameter algorithm is 
usually applied when the parameter is small. In this situation the exponential compo- 
nent of the runtime becomes a relatively small multiplicative or additive constant, that 
is the problem is solved by a polynomial (usually even a low polynomial) algorithm. 
Nevertheless, even if the parameter is small, the researchers often prefer to use impre- 
cise heuristic methods simply because they are faster. We argue that in this situation 
a fixed-parameter algorithm may be still of a considerable use because it can help to 
evaluate the quality of the considered heuristics. In particular, in our experiments the 
use of a fixed-parameter algorithm allowed to observe that the heuristic SGA80(DFS) 
almost always returns an optimal solution. That is, although the heuristic is the slowest 
among the considered ones, this is still not the reason to discard it: the heuristic can be 
the best choice when the quality is particularly crucial. Thus we represent a novel way 
of application of fixed-parameter algorithms which, we believe, would be a significant 
contribution to the applied research related to fixed-parameter computation. 

To design a fixed-parameter algorithm for the DMERN problem we in fact design 
a fixed-parameter algorithm for the maximum balanced subgraph problem using the 
equivalence result from ifPD . The fixed-parameter algorithm for the latter problem 
is designed by showing its equivalence to the bipartization problem and then using a 
fixed-parameter algorithm for the bipartization problem first proposed in ll25l and then 
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refined in |fl9l . 

The rest of the paper is organized as follows. In Section|2]we introduce necessary 
notation, Section[3]presents the SGA heuristic and its variants, and Section|4]introduces 
the fixed-parameter algorithms. Section|5]describes a fixed-parameter algorithm for the 
maximum balanced subgraph problem. In Section [6] we report empirical results and 
analyze them. Concluding remarks are made in Section|7] 

2 Embedded networks and signed graphs 

In this section, we assume, for simplicity, that A is a (0, ±l)-matrix itself (since all rows 
containing entries not from the set { — 1, 0, +1} cannot be part of a reflected network). 
Here we allow graphs to have parallel edges, but no loops. A graph G — (V, E) along 
with a function s : E— >{ — , +} is called a signed graph. Signed graphs have been 
studied by many researchers, see, e.g., Ifl5l [161 [T71 l26l . 

We assume that signed graphs have no parallel edges of the same sign, but may 
have parallel edges of opposite signs. An edge is positive (negative) if it is assigned 
plus (minus). For a (0, ±l)-matrix A — [aik] with n rows, we construct a signed 
graph G(A) as follows: the vertex set of G(A) is {1, 2, . . . . n}; G(A) has a positive 
(negative) edge ij if and only if aik = —a>jk ^ (aik = ajk ^ 0) for some k. Let 
G = (V, E, s) be a signed graph. For a non-empty subset W of V, the W -switch of G 
is the signed graph G w obtained from G by changing the signs of the edges between 
W and V (G) \ W . A signed graph G = (V, E, s) is balanced if there exists a subset 
W of V (W may coincide with V) such that G w has no negative edges. Let 77(G) be 
the largest order of a balanced induced subgraph of G. 

The following important result was proved in ff3l . This result allows us to search 
for a largest balanced induced subgraph of G(A) instead of a largest reflected network 
in A. 

Theorem 1. H13V Let A be a (0, ±l)-matrix. A set R of rows in A forms a reflected 
network if and only if the vertices of G(A) corresponding to R induce a balanced 
subgraph of G (A). In particular, v(A) = rj(G(A)). 

3 SGA and its Variations 

The heuristic SGA introduced in |fl3ll is based on the following: 
Lemma 1. H13V Every signed tree T is a balanced graph. 

Proof. We prove the lemma by induction on the number of edges in T. The lemma 
is true when the number of edges is one. Let a; be a vertex of T of degree one. By 
the induction hypothesis, there is a set W C V(T) — x such that (T — x) w has no 
negative edges. In T w the edge e incident to x is positive or negative. In the first case, 
let W = W and the second case, let W — W U {x}. Then, T w has no negative 
edges. □ 
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Figure 1: Illustration for SGA; M is the subgraph G^ induced by the negative edges 
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Heuristic SGA: 

Step 1: Construct signed graph G = G(A) = (V, E, s). 
Step 2: Find a spanning forest T in G. 

Step 3: Using a recursive algorithm based on the proof of LemmaQ] compute 
W C V such that T w has no negative edges. 

Step 4: Let N be the subgraph of G w induced by the negative edges. Apply 
the following greedy-degree algorithm [23] to find a maximal independent set 
/ in N: starting from empty /, append to / a vertex of N of minimum degree, 
delete this vertex together with its neighbors from N, and repeat the above 
procedure till N has no vertex. 
Step 5: Output /. 



For a graph H with connectivity components Hi, ... , H p , a spanning forest is the 
union of spanning trees Ti, . . . , T p of Hi , . . . , H p , respectively. The second step of the 
algorithm (i.e., finding of spanning forest) can be implemented in a number of different 
ways. We tried the following ones: 

1. Random Search (RS) starts from marking some vertex. On every iteration, it 
finds some edge (positive or negative; double edges are not considered) between 
a marked vertex v and an unmarked vertex u. It marks u and includes or excludes 
it from the set W according to the sign of the edge uv. If no edges between 
marked and unmarked vertices are found, the algorithm marks some random 
vertex again. This was the method of computing a spanning forest used for the 
experiments reported in JT3). 

2. BFS is a well known algorithm for constructing of spanning forests. It maintains 
a FIFO queue. On every iteration, it takes a vertex from this queue, marks all its 
unmarked neighbors into the queue and adds these neighbors to the queue. If the 
queue is empty, some unmarked vertex is marked and added into the queue. In 
our implementation, we take a vertex of maximum degree in this case. 

3. DFS is another well known algorithm for constructing of spanning forests. It 
starts from some vertex and then calls itself recursively for every of its neighbors 
which still are not included in the spanning forest. 



Proposition 1. [13]IfG is balanced, then I = V. 
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Proof. It is well-known (see, e.g., Theorem 2.8 in [13]) that a signed graph is balanced 
if and only if it does not contain cycles with odd number of negative edges. Let T be 
a a spanning forest in G. Since T w has no negative edges, G w cannot have negative 
edges. Indeed, if xy was a negative edge in G w , it would be the unique negative edge 
in a cycle formed by xy and the (x, y)-path of T w , a contradiction. □ 

Gutin and Zverovitch lfl4l investigated a repetition version of SGA where Steps 
2-4 were repeated several time (each time the vertices of G were pseudo-randomly 
permuted and a new spanning forest of G was built). They found out that three times 
repetition of SGA gives about 1% improvement, while 80 times repetition of SGA 
leads to 2% improvement, on average. In our experiments we used a larger test bed 
and better scaling procedure than in [ 14 1 and, thus, we run SGA and its 3 and 80 times 
repetitions on the new set of instances of the DMERN problem (see Section [6j. We 
will denote these repetition versions of SGA by SGA3 and SGA80, respectively. 

In Section [6] we also report results on another modification of SGA, SGA+VC, 
where we replace Step 4 with finding a vertex cover C of G w and setting / = V(G W )\ 
C. Since the vertex cover problem is well studied in the area of parameterized complex- 
ity lfTl[8l l22l . to find C we can use a fixed-parameter algorithm for the problem. 

4 Fixed-Parameter Algorithmics 

We recall some most basic notions of fixed-parameter algorithmics (FPA) here, for a 
more in-depth treatment of the topic we refer the reader to the monographs [10 TT1I221 . 

FPA is a relatively new approach for dealing with intractable computational prob- 
lems. In the framework of FPA we introduce a parameter fc, which is often a positive 
integer (but may be a vector, graph, or any other object for some problems) such that the 
problem at hand can be solved in time 0(f(k)n c ), where n is the size of the problem 
instance, c is a constant not dependent on n or k, and f(k) is an arbitrary computable 
function not dependent on n. The ultimate goal is to obtain /(fc) and c such that for 
small or even moderate values of k the problem under consideration can be completely 
solved in a reasonable amount of time. 

As an example, consider the Vertex Cover problem (VC): given an undirected graph 
G (with n vertices and m edges), find a minimum number of vertices such that every 
edge is incident to at least one of these vertices. In the (naturally) parameterized version 
of VC, fc-VC, given a graph G, we are to check whether G has a vertex cover with at 
most k vertices. fc-VC admits an algorithm of running time 0(1.2738 fc + kn) obtained 
in [8 1 that allows us to solve VC with k up to several hundreds. Without using FPA, we 
would be likely to end up with the obvious algorithm of complexity 0(mn k ). The last 
algorithm is far too slow even for small values of fc such as fc = 10. 

Parameterized problems that admit algorithms of complexity 0(f(k)n c ) (we refer 
to such algorithms as, fixed-parameter) are called fixed-parameter tractable (FPT). No- 
tice that not every parameterized problem is FPT, but there are many problems that are 

fpt nnummi. 
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5 Minimum Balanced Deletion problem 

By our discussions above, we are interested in the following parameterized problem. 



The minimum balanced deletion problem (MBD) 

Input: A signed graph G = (V, E, s), an integer k. 
Parameter: k. 

Output: A set of at most k vertices whose removal makes G balanced or 'NO' 
if no such set exists. 

We show that the MBD problem is FPT by transforming it into the Bipartization 
problem defined as follows. 

The Bipartization problem 

Input: A graph G, an integer k 
Parameter: k 

Output: A set of at most k vertices whose removal makes G bipartite or 'NO' 
if no such set exists. 

The transformation is described in the following theorem. 

Theorem 2. The MBD problem is FPT and can be solved in time O* (3 fe ). 

Proof. It is well-known (see, e.g., Theorem 2.8 in [13]) that a signed graph is balanced 
if and only if it does not contain cycles involving odd number of negative edges. Hence, 
the MBD problem in fact asks for at most k vertices whose removal breaks all cycles 
containing an odd number of negative edges. 

Let G' be the (unsigned) graph obtained from G by subdividing each positive edge. 
In other words, for each positive edge {u, v}, we introduce a new vertex w and replace 
{u,v}by {u,w} and {w, v}. We claim that G has a set of at most k vertices breaking 
all cycles with an odd number of negative edges if and only if G' can be made bipartite 
by removal of at most k vertices. 

Assume the former and let K be a set of at most k vertices whose removal breaks all 
cycles with an odd number of negative edges. It follows that G' — K is bipartite. Indeed, 
each cycle C' of G' — K can be obtained from a cycle C of G — K by subdivision of 
its positive edges. Hence, C' can be of an odd length only if C has an odd number of 
negative edges which is impossible according to our assumption about K. 

Conversely, let if be a set of at most k vertices such that G' — K is bipartite. We 
may safely assume that K does not contain the new vertices subdividing positive edges: 
otherwise each such vertex can be replaced by one of its neighbors. Thus, K C V(G). 
Observe that G — K does not have cycles with odd number of negative edges. Indeed, 
by subdividing positive edges, any such cycle translates into an odd cycle of G 1 — K 
in contradiction to our assumption about K. 

It follows from the above argumentation that the MBD problem can be solved as 
follows. Transform G into G 1 and run on G 1 the 0*(3 fc ) algorithm solving the bi- 
partization problem |fl9l . If the algorithm returns 'NO' then return 'NO'. Otherwise, 
replace each subdividing vertex by one of its neighbors and return the resulting set of 
vertices. Clearly, the complexity of the resulting algorithm is O* (3 fc ). □ 
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Remarks. The usual trick to avoid the subdivided vertices to be selected to the 
resulting solution would be to make k + 1 copies of each vertex. However, such ap- 
proach would increase the runtime of the resulting implementation and hence we have 
used a slightly more sophisticated method. Unfortunately, it is not known yet whether 
the Bipartization problem has a polynomial-size problem kernel |19j. Thus, it is not 
known yet whether the MBD problem has a polynomial-size problem kernel. (If one 
was known, we could try to use it to speed up our fixed-parameter algorithm.) 

Note that a version of the MBD problem, where edge-deletions rather than vertex- 
deletions are used was considered in ||20l l9l. 

6 Experimental Evaluation 

In this section we provide and discuss our experiment results for the heuristics SGA, 
SGA3, SGA80, SGA+VC descried in Section[3]and the exact algorithm given in Sec- 
tion [5] Note that in our experiments we use a larger test bed and better scaling proce- 
dure than in lfT4l . 

6.1 Scaling Procedure 

Recall that we consider an LP problem in the standard form stated as 

Minimize {p T x; subject to Ax — b, x > 0}. 

In Section|2] to simplify our notation we assumed that A is a (0, ±l)-matrix. However, 
in general, in real LP problems A is not a (0, ±l)-matrix. Therefore, in reality, the first 
phase in solving the DMERN problem is applying a scaling procedure whose aim is to 
increase the number of (0, ±l)-rows by scaling rows and columns. Here we describe 
a scaling procedure that we have used. Our computational experiments indicate that 
this scaling is often better than the scaling procedures we found in the literature. Let 
us describe our scaling procedure. Let A — [aij] nxm . 

First we apply simple row scaling, i.e., scale all the rows which contain only zeros 
and ±2, where x > is some constant: for every i £ {1, 2, . . . , n} set ay = aij/x for 
j = 1, 2, . . . ,mif ay £ {0, —x, +x} for every j e {1, 2, ... , m}. 

Then we apply a more sophisticated procedure. Let [ri] n be an array of boolean 
values, where ri indicates whether the ith row is a (0, ±l)-row. Let [bj) m be an array 
of boolean values, where bj indicates whether the jth column is bounded, i.e., whether 
it has at least one nonzero value in a (0, ±l)-row: for some j e {1,2,..., m} the value 
bj = true if and only if there exists some i such that ri — true and ay ^ 0. 

Next we do the following for every non (0, ±l)-row (note that at this stage any 
non (0, ±l)-row contains at least two nonzero elements). Let J be the set of indices 
of bounded columns with nonzero elements in the current row c: J — {j : a C j ^ 
and bj — true}. If J = 0, i.e., all the columns corresponding to nonzero elements 
in the current row are unbounded, then we simply scale every of these columns: ay — 
dij/a C j for every i = 1,2, ... ,n and for every j such that a C j ^ 0. If J ^ and 
a C j £ {+x, —x} for every j G J, where x is some constant, then we scale accordingly 
the current row (a C j — a C j/x for every j S {1,2,..., to}) and scale the unbounded 
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columns: a«j = aij/a c j for every j : ^ J if a c j ^ 0. Otherwise we do nothing for the 
current row. 

Every time when we scale rows or columns we update the arrays r and b. 

Since the matrices processed by this heuristic are usually sparse, we use a special 
data structure to store them. In particular, we store only nonzero elements providing 
the row and column indices for each of them. We also store a list of references to the 
corresponding nonzero elements for every row and for every column of the matrix. 

6.2 Computational Experience 

The computational results for all heuristics apart from SGA+VC as well as for the ex- 
act algorithm are provided in TableQ] As a test bed we use all the instances provided in 
Netlib ( http : / /net lib . org/lp/data/) . All algorithms were implemented in 
C++ and the evaluation platform is based on an AMD Athlon 64 X2 3.0 GHz processor. 

For the exact algorithm we used a code of Huffnei fFittp : //theinf 1 . inf ormatik . uni- jena . de/ ~huef f ner 
In SGA+VC we used a vertex cover code based on |2|. Since the exact algorithm can 
potentially take a very long time, we introduced a timeout: if after one hour of work 
the algorithm is unable to compute the output, it terminates. 

Denote by n the number of (0, ±l)-rows in the instance, i.e., the number of vertices 
in the corresponding signed graph G. The column k reports the difference between n 
and the number of vertices in a maximum induced balanced subgraph of G found by the 
exact algorithm. The rows where k is not given correspond to the instances where the 
algorithm terminated after one hour without computing the output. The nine columns 
following column k report the same differences found by heuristics SGA, SGA3, and 
SGA80. The first three columns are related to SGA. The columns RS, BFS, and DFS 
are related to the way of computing the spanning tree (see Section 3 for the detailed 
explanation). The next three columns are related to SGA3 and the last three columns 
are related to SGA80 with the analogous meaning of particular columns. The column 
t reports the runtime taken by the exact algorithm. The time provided in the column 
V is the average time required for SGA(RS), SGA(BFS) and SGA(DFS)Q to proceed 
once. Our experiments show that for SGA3 this time is approximately 3 times larger 
and for SGA80 80 times larger. It also appears that there is no significant difference 
between running times of SGA based on RS, BFS or DFS. 

The 'Average' row shows the average value of the respective columns. The Avg. 
diff.' row shows how far, on average, is a particular modification of SGA from the 
optimal solutions. '# exact sol.' shows the number of instances for which a particular 
modification of SGA obtained an optimal solution. Both Avg. diff.' and '# exact sol.' 
consider only the instances with the knows optimal solutions. 



the acronym in the parenthesis correspond to the algorithm of computing the spanning tree. 
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Table 1: Experiment results. 



Instance 


K 




SGA 




SGA3 




SGA80 


Time, s 


RS 


BFS 


DFS 


RS 


BFS 


DFS 


RS BFS 


DFS 


t 1 1 


25FV47 


15 


25 


TO 

38 




25 


38 


1 

18 


22 


38 


16 


4.40 0.004 


80BAU3B 


— 


42 


ACt 




40 


AO 


At\ 


40 


40 


Ar\ 


> lh 0.106 


ADLITTLE 


1 


1 


1 


l 


1 


1 


1 


1 


1 


1 


0.02 0.000 


AFIRO 
































0.00 0.000 


AGG 


— 


107 


108 


104 


106 


108 


104 


104 


108 


103 


> lh 0.001 


AGG2 


— 


85 


91 


83 


85 


91 


83 


83 


91 


83 


> lh 0.001 


AGG3 


— 


85 


91 


83 


85 


91 


83 


83 


91 


83 


> lh 0.001 


BANDM 


23 


24 


24 


24 


23 


24 


23 


23 


24 


23 


1493.12 0.001 


BEACONFD 


3 


3 


3 


3 


3 


3 


3 


3 


3 


3 


0.00 0.000 


BLEND 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


0.00 0.000 


BNL1 


14 


17 


19 


15 


17 


18 


14 


14 


18 


14 


1.83 0.003 


BNL2 


— 


127 


109 


103 


110 


107 


101 


99 


105 


86 


> lh 0.061 


BOEING 1 


— 


49 


61 


42 


49 


61 


42 


48 


60 


42 


> lh 0.001 


BOEING2 


15 


17 


17 


16 


17 


17 


16 


15 


17 


15 


0.05 0.000 


BORE3D 


12 


14 


15 


13 


14 


15 


12 


12 


15 


12 


0.14 0.001 


BRANDY 


6 


7 


8 


6 


6 


8 


6 


6 


8 


6 


0.00 0.000 


CAPRI 


— 


40 


43 


37 


37 


40 


34 


34 


40 


33 


> lh 0.001 


CYCLE 


— 


34 


36 


34 


34 


36 


34 


34 


36 


34 


> lh 0.020 


CZPROB 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


0.27 0.018 


D2Q06C 


— 


67 


97 


67 


67 


96 


67 


67 


94 


65 


> lh 0.042 


D6CUBE 


— 


61 


50 


50 


52 


50 


50 


46 


47 


50 


> lh 0.003 


DEGEN2 


— 


234 


237 


219 


234 


237 


219 


226 


234 


218 


> lh 0.011 


DEGEN3 


— 


822 


819 


769 


822 


819 


769 


815 


819 


769 


> lh 0.197 


DFL001 


— 


2818 2997 2603 


2818 2956 2603 


2802 


2903 2585 


> lh 1.871 


E226 


15 


18 


19 


15 


17 


19 


15 


16 


19 


15 


1.09 0.001 


ETAMACRO 


12 


20 


31 


19 


20 


31 


19 


20 


26 


16 


0.47 0.001 


FFFFF800 


— 


50 


69 


46 


50 


69 


46 


41 


65 


39 


> lh 0.002 


FINNIS 


— 


121 


127 


120 


121 


127 


120 


119 


127 


119 


> lh 0.003 


FIT ID 


6 


6 


7 


6 


6 


7 


6 


6 


7 


6 


0.00 0.000 


FIT1P 
































0.00 0.000 


FIT2D 


6 


7 


7 


6 


6 


7 


6 


6 


7 


6 


0.00 0.004 


FIT2P 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


0.00 0.007 


FORPLAN 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


0.00 0.000 


GANGES 




83 


84 


82 


83 


84 


82 


77 


84 


77 


> lh 0.021 


GFRD-PNC 




68 


68 


95 


68 


68 


86 


68 


68 


80 


> lh 0.008 


GREENBEA 




48 


60 


46 


48 


60 


46 


45 


57 


41 


> lh 0.043 


GREENBEB 




48 


60 


46 


48 


60 


46 


45 


57 


41 


> lh 0.044 


GROW 15 
































0.00 0.000 


GROW22 
































0.00 0.000 


GROW7 
































0.00 0.000 


ISRAEL 


8 


9 


10 


8 


8 


10 


8 


8 


10 


8 


0.02 0.000 


KB 2 


1 


1 


3 


2 


1 


3 


1 


1 


3 


1 


0.00 0.000 


LOTFI 


18 


24 


24 


19 


22 


20 


19 


19 


20 


19 


11.23 0.001 


MAROS-R7 
































0.00 0.011 


MAROS 


11 


17 


14 


21 


15 


14 


18 


12 


14 


11 


0.23 0.005 


MODSZK1 




237 


267 


237 


237 


267 


237 


237 


267 


235 


> lh 0.004 


NESM 


10 


13 


13 


13 


11 


13 


11 


10 


13 


10 


0.03 0.003 


PEROLD 




28 


29 


25 


26 


29 


24 


24 


27 


23 


> lh 0.003 


PILOTJA 


16 


18 


19 


18 


17 


19 


16 


16 


19 


16 


11.72 0.004 


PILOT 




45 


45 


44 


42 


44 


41 


41 


44 


41 


> lh 0.008 
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SGA SGA3 SGA80 Time, s 



lll,>lullLU 


k 


RS 


BFS 


DFS 


RS BFS 


DFS 


RS 


BFS 


DFS 


t t\ 


PILOT.WE 


— 


34 


36 


31 


30 


36 


30 


28 


36 


27 


> lh 0.003 


PILOT4 


3 


3 


3 


3 


3 


3 


3 


3 


3 


3 


0.00 0.001 


PILOT87 


— 


77 


87 


74 


76 


87 


71 


70 


87 


69 


> lh 0.015 


PILOTNOV 


19 


21 


24 


21 


21 


24 


20 


19 


24 


19 


201.29 0.005 


RECIPE 
































0.00 0.000 


SC105 


16 


17 


32 


22 


17 


31 


21 


17 


31 


19 


12.56 0.000 


SC205 


— 


36 


68 


50 


36 


67 


47 


36 


67 


41 


> lh 0.000 


SC50A 


8 


8 


12 


8 


8 


11 


8 


8 


11 


8 


0.02 0.000 


SC50B 


6 


6 


9 


8 


6 


9 


6 


6 


9 


6 


0.02 0.000 


SCAGR25 
































0.03 0.002 


SCAGR7 
































0.02 0.000 


SCFXM1 


12 


13 


14 


13 


12 


14 


13 


12 


14 


12 


0.30 0.001 


SCFXM2 


— 


26 


28 


26 


26 


28 


26 


24 


28 


26 


> lh 0.003 


SCFXM3 


— 


39 


42 


39 


38 


42 


39 


36 


42 


39 


> lh 0.006 


SCORPION 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


0.02 0.001 


SCRS8 


9 


9 


9 


9 


9 


9 


9 


9 


9 


9 


0.06 0.002 


SCSD1 
































0.03 0.000 


SCSD6 
































0.14 0.001 


SCSD8 
































0.59 0.001 


SCTAP1 
































0.00 0.000 


SCTAP2 
































0.00 0.006 


SCTAP3 
































0.00 0.012 


SEBA 


— 


274 


285 


267 


274 


285 


267 


269 


285 


265 


> lh 0.021 


SHARE IB 


4 


5 


4 


4 


4 


4 


4 


4 


4 


4 


0.00 0.000 


SHARE2B 


6 


6 


6 


6 


6 


6 


6 


6 


6 


6 


0.00 0.000 


SHELL 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


1.51 0.006 


SHIP04L 


— 


36 


36 


36 


36 


36 


36 


36 


36 


36 


> lh 0.006 


SHIP04S 


— 


36 


36 


36 


36 


36 


36 


36 


36 


36 


> lh 0.005 


SHIP08L 




64 


64 


64 


64 


64 


64 


64 


64 


64 


> lh 0.023 


SHIP08S 


— 


64 


64 


64 


64 


64 


64 


64 


64 


64 


> lh 0.015 


SHIP12L 




96 


96 


96 


96 


96 


96 


96 


96 


96 


> lh 0.046 






96 


96 


96 


96 


96 


96 


96 


96 


96 


v lh (1 ()-2() 
^> 111 \J.\J3\J 


SIERRA 




400 


455 


420 


400 


455 


397 


387 


444 


388 


> lh 0.030 


STAIR 


8 


11 


12 


9 


9 


12 


9 


8 


12 


8 


0.02 0.000 


STANDATA 




53 


57 


59 


53 


57 


59 


53 


57 


59 


> lh 0.002 


STANDGUB 




53 


57 


59 


53 


57 


59 


53 


57 


59 


> lh 0.002 


STANDMPS 




54 


58 


73 


54 


58 


73 


54 


58 


70 


> lh 0.004 


STOCF0R1 
































0.00 0.000 


STOCFOR2 




258 


311 


202 


258 


307 


202 


243 


305 


197 


> lh 0.050 


TUFF 


16 


26 


16 


18 


17 


16 


16 


16 


16 


16 


0.58 0.001 


VTP.BASE 


4 


6 


6 


6 


4 


6 


4 


4 


6 


4 


0.00 0.000 


WOOD IP 
































0.02 0.001 


WOODW 
































0.00 0.007 


Average 
Avg. diff. 
# exact sol. 




79.3 84.85 75.57 
1.28 2.15 0.93 

33 31 36 


78.55 
0.78 
40 


84.2 74.82 
2.02 0.52 
31 44 


76.91 83.19 73.54 
0.35 1.93 0.17 
48 31 50 


32.26 0.030 



Observe that the exact algorithm completed its computations for 54 instances out 
of the total of 93, and for 52 instances the running time was at most 1 minute. Note 
that SGA(DFS) achieved the optimal solution in 36 out of 54 cases, SGA3(DFS) in 44 
cases and SGA80(DFS) in 50 cases. 

The DFS-based version of SGA clearly outperforms the RS- and BFS-based ver- 
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sions. This is also true for SGA3 and SGA80. We are not able to justify this result but 
it is strongly confirmed by our extensive experimentation. 

The results with SGA(DFS)+VC are not provided in Table[TJsince SGA(DFS)+VC 
managed to improve SGA(DFS) only for 4 instances: 



Instance 


k 


fcsGA(DFS)+VC 


fcsGA(DFS) 


fcsGA80(DFS) 


BOEING2 


15 


15 


16 


15 


DEGEN2 




218 


219 


218 


DEGEN3 




764 


769 


769 


DFL001 




2599 


2603 


2585 



SGA(DFS)+VC was not able to proceed in 30 minutes for AGG, AGG2, AGG3, 
FINNIS, MODSZK1 and SIERRA. Note that /jsga(dfs)+vc < &sga80(dfs) for onl y one 
instance while fcsGA(DFS)+vc > ^sgaso(dfs) for 39 instances and recall that for 6 in- 
stances SGA(DFS)+VC did not terminate in the given time. Since the running time of 
SGA(DFS)+VC usually exceeds that of SGA80(DFS) and the quality of SGA(DFS)+VC 
is not much different even from that of SGA(DFS), SGA+VC appears to be of little 
practical interest. However, SGA+VC demonstrates that there is no need to replace 
Step 4 of SGA by a more powerful heuristic or exact algorithm and that in order to 
improve results one should pay attention to the Step 2 of SGA. In other words, this 
negative result pointed out the proper area for further research. Indeed, replacement 
of RS algorithm for the spanning tree building with DFS dramatically improved the 
heuristic quality. 

Although SGA80 is slower than SGA and SGA3, it is much more precise and 
its running time is still reasonable. It follows that SGA80(DFS) is the best choice 
with respect to the tradeoff between running time and solution quality. Observe that in 
almost all (50 out of 54) the cases feasible for the exact algorithm, SGA80(DFS), being 
much faster than the exact algorithm, managed to compute an optimal solution! Note 
that this conclusion can be made only given the knowledge about the optimal solution 
and without the considered fixed-parameter algorithm such knowledge would be very 
hard to obtain (for example, the instance PILOTNOV with n = 329 and fc m j n = 19 
would hardly be feasible to a brute-force exploration of all (^g 9 ) possibilities). 

7 Conclusions 

One of contributions of this paper is demonstrating a novel way of use of fixed-parameter 
algorithms where they do not substitute heuristic methods but are used to evaluate them. 
As a case study, we considered heuristics for the problem of extracting a maximum- 
size reflected network in an LP problem. The use of fixed-parameter algorithms helped 
us to investigate a state-of-the-art heuristic, improve it and allowed us to arrive at an 
interesting observation that the improved version of the considered heuristic almost 
always returns an optimal solution. 
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We believe that this way of applying fixed-parameter algorithms can be useful for 
other problems as well. One candidate might be the problem of finding whether the 
given CNF formula has at most k variables so that their removal makes the resulting 
formula Renameable Horn. This is called the Renameable Horn deletion backdoor 
problem and was recently shown FPT [24]. Heuristics for this problem are widely 
used in modern SAT solvers for identifying a small subset of variables on which an 
exponential-time branching is to be performed lETI . Currently it is unclear whether 
substituting a heuristic approach by the exact fixed-parameter algorithm would result 
in a better SAT solver. But even if it is not the case, the exact algorithm can be still of 
a considerable use for ranking the heuristic techniques, especially as producing small 
Renameable Horn backdoors is vitally important for reducing the exponential-time im- 
pact on the runtime of SAT solvers. 
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